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1. INTRODUCTION 


T ^ rUSt VeCt0C1 " 8 h3S *- -***- - - — * t. enhance maneuver- 

ablUty and Perf °— • tact i.al aircraft. This concept 

usually Involves the installation o£ a „ ltlfunction no2Ue _ mch the ^ 

tha engine thrust can be changed in faction without changing the attitude ol 

tha alrCraft ' *“ “ - thrust induces a significant change 

in the aerodynamic forces acting on the aircraft rh p 

6 aircraft. Therefore, this device can 

lift-augmenting as wen as stabllity and controi purpogeg> ^ 

'"a tHrUSt 13 d6fleCted ^ the ^tudinal direction, the lift force and 
Pitching stahility can he manipulated, while yawing stability can be 
controlled by directing the thrust in the lateral direction. 

Experimental inyestigation of the thrust vectoring concept has been under 
»ay for decades (Refs. 1-3). However, computational methods have not Rept 
Pure with experiments. One computational method based on a panel method was 

loped in Ref. 4. For configurations having multifunction nozzles, the 

no2zle ex . t becQmes compUcated> jet entrainment ^ 

well as shock wave induced by the jet stream will be preS ent Also 

present. Also, pressure 

difference across the let exit- vtn 

j exit will produce vortices. Due to the limitation 

of flow singularity methods, these critical flow characteristics cannot be 

predicted by a panel code. Another drawback of panel methods is the absence 

of. energy and temperature terms -fn m, d 

terms in the computation. Also, jet exit conditions 

cannot be included in the calculation. Therefore the i t 

neretore, the jet stream is frequent- 

ly modeled as a flap or wahe. To fully understand the induced aerodynamic 

affects from the jet stream deflection, methods other than a panel approach 
would be needed. 


1 


Methods in computational fluid dynamics (CFD) are more suitable for the 
present problem. CFD applications can be described as applying the laws of 
flow physics to mathematical equations and solving them with a computer. The 
concept of CFD can be traced back to the invention of digital computers. 
However, applications to complicated problems did not succeed until 
recently. In this decade, successful methods have been developed in scheme 
stability, solution techniques, grid generation, etc. Thanks to the progress 
of CFD techniques, flow properties related to thrust vectoring such as shock 
wave and vorticity can be captured by various methods. Also, with the 
development of numerical grid generation, complicated configurations can be 
modeled accurately. 

Two CFD methods, coupled with a numerical grid generation program, have 
been used in this study. These two CFD methods, Euler and thin-layer Navier- 
Stokes algorithms (Refs. 5, 6), both utilize the finite volume scheme, which 
is suitable to model the jet momentum transfer. The numerical grid generation 
program EAGLE (Ref. 7), with the options of algebraic and elliptic grid 
generation systems, has been used to produce grids in this study. 

Theoretical approaches of these three programs will be summarized in this 
report. Some computational results will also be presented. 


2 



2. THEORETICAL APPROACH 


Generally, tasks in developing CFD programs include derivation of governing 
equations, solution of the flow equations, accuracy and stability analysis of 
the scheme, grid generation, etc. In this section, theoretical outlines and 
comparisons of the two programs, FLO-59 and CFL3D, will be discussed. 

The program FLO-59 solves Euler equations, which feature conservation of 
mass, momentum, and energy In inviscid flow. Because of the omission of 
viscous terms, this program has the advantage of saving computer resources and 
CPU time. However, since this formulation does not Include viscous effects, 
phenomena such as flow separation and secondary vortex cannot be predicted 
accurately. The program CFL3D utilizes thin-layer Navier-Stokes equations, 
thus more detailed flow properties can be obtained. Yet this program requires 
more CPU time and memory space. From experience, choice of these two programs 
depends on the purpose of applications. For example, Euler codes have great 
success In predicting some primary vortex features. It is mainly due to the 
fact that the primary vortex Is generated by pressure differences between two 
regions. In predicting secondary vortex flow, Navier-Stokes equations are 
needed because the secondary flow Is caused by viscous effect. Therefore, if 
the flow phenomena to be predicted are caused by viscous effect, Navier-Stokes 
equations should be used. Otherwise, the Euler approach should be able to 
provide a satisfactory solution. 

In terms of solving governing equations, both programs use the finite 
volume approach. In both cases, the accuracy of the solution depends on the 
smoothness of the grids. 

In discretizing a differentiation into finite difference forms, central 
difference Is usually used. This formulation ensures a second order accu- 
racy. From the consideration of accuracy, a central difference Is better 
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than either the backward or forward difference. However, in some particular 
situation, backward difference has to be used to satisfy the physics of the 
flow field. For example, in a supersonic flow region, the Mach lines drawn 
from a point divide the flow into zones of influence. Properties at this 
point can only be affected by flow in the zone of influence. Therefore, while 
formulating differentiation with difference forms, only the points in the zone 
of influence should be taken into consideration. In this situation, the 
backward difference has to be used. For subsonic flow, characteristics of 
each point can be affected by any other points in the flow field. Thus, 
central difference is suitable for subsonic flow calculation. 

In transonic flow computation, the formulation of the differential equa- 
tions becomes more difficult due to the combination of sub- and supersonic 
regions. As stated earlier, in the subsonic region, the central difference 
should be used. On the other hand, In the supersonic region, backward 
difference has to be used. Upwind difference was developed to solve this 
problem. Take the differentiation dF/d£ as an example. The formulation can 
be written as 


dF 

as 


F - F F - F 

, i + 1 i-U . ,, w i i-U 

777 ) + (1 “ H>< 7? ) 


2AS 


AS 


\i =* 1 subsonic (1) 

\i = 0 supersonic 

Therefore, central difference is used in the subsonic region. In a supersonic 
region, backward difference is used. 

FLO-59 adds dissipative terms to the Euler equations to achieve upwind 
differencing. On the other hand, CFL3D has the options of flux difference 
splitting and flux vector splitting to separate the flow field Into different 
regions. Details of upwind differencing of these two programs will be 
discussed later. 
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The most distinct difference between these two programs is the single- 
block and multiblock operations. FLO-59 utilized the single-block concept. 
Therefore, applications are limited to simple configurations. On the other 
hand, CFL3D has the capability of multiblock computation. The development of 
multiblock computation enables CFD applications to complicated configurations. 
Also, this concept makes the grid generation an easier task. For example, an 
individual block can be formed of the region where smooth continuation of the 
grids is not possible. Therefore, smooth grids can be constructed in each 
block. Flow properties can also be solved accurately in each block. Trans- 
mission of flow characteristics among blocks can be achieved through interpo- 
lations. Even the grid spacings are different among blocks. 

2.1 Euler Method (FLO-59) 

2.1.1 Governing Equations 

The equations can be expressed in an integral form as follows. 

f- /// q dQ + JJ • n dS = 0 (2) 

<3t Q S 

where 

]5 ~ N e +Ne +Ne 
xx y y z z 




f pu N 
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f \ 

pw 
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puv 


pwu 
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> N x - 

puv 

* N y 

pv 2 +p 

, N z - 

pwv 

pw 
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pvw 
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2.1*2 Finite Volume Sche: 


These Euler equations can be discretized using the finite volume 
concept. The physical domain is divided into hexahedral cells. At each 
finite volume cell, the Euler equations can be written as 


dt ^ijk^ijk^ + ^ijk^ijk ^ 


( 4 ) 


where is the cell volume and is the surface integral. The second 

term ± j 1c re P resents the net flux out of the cell which is balanced by the 

rate of change of q in the cell volume. The net flux term is given by 


F ijk q ijk = N i+l/2,j,k S i+l/2,j,k “ N i-l/2,j,k S i-l/2,j,k 

+ N i,j+l/2,k S i,j+l/2,k " N i,j-l/2,k S i,j-l/2,k 

+ N i,j,k+l/2 S i,j,k+l/2 " N i,j,k-l/2 S I,j,k-l/2 
where S denotes the area of the cell surface. 


2.1.3 Upwind Differencing 


The use of central difference in Equation (5) ensures that this scheme is 
second order accurate in case the grid spacing is uniform. However, this 
scheme is suitable for subsonic flow computation only. For supersonic flow, 
backward differencing should be used in the discretization. An alternative 
form of Equation (4) can be written as 


d__ 

dt 


^ J ijk q ijk^ + F ijk q ijk 


°ijk q ijk 


0 


( 6 ) 


where 


Dq » D q + D q, + D q 
n ijk xMjk y ijk z^ijk 

D x q ijk " d i+l/2,j,k " d i-l/2,j,k 

°y q ijk d i,j+l/2,k d i,j-l/2,k 
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( 8 ) 


°z q ijk " d i,j,k+l/2 " d i,j,k-l/2 


i+l/2,i,k 


111LL1A [e <2) ( 


At 


i + l/2,j ,k M i+l,j ,k H i,j,k 


) 


- £ 


i+l/2,j,k (q i+2,j,k' 3q i+l,j,k + 3q ijk " q i+j,k 


)] 


( 9 ) 


P - 2 P + p 

define v... = _±±L 1 A _LiiS _i±LA. 

Ip + 2 P + p I 

1 zt ijk + p i+j,k' 


( 2 ) 

C l+l/2,j,k- k -“<Vl,J,k> V 5 


.( 2 ) 

i+1, 

- max[0 (k^ - 31 

e I+l/2,j,k maxlU,U e i+l/2,j,k )1 


( 10 ) 


( 11 ) 


( 12 ) 


In Equation ( 11 ), the artificial sensor v, .. is used to sense the existence of 

lJK 

shock wave. Theoretically, the pressure distribution is smooth in the flow 
field if no shock wave exists. Therefore, v , is close to zero. Central 
difference is used in this situation. If the shock occurs In the flow field, 
the pressure rises across the shock. The sensor is not zero; thus, backward 
difference is brought into the computation. 


2.1.4 Time— Stepping Scheme 


At each point. Equation ( 5 ) is written as 

+ J (Fq - Dq) = 0 (13) 

For the integration of Equation (6), the time~stepping scheme used is a 
fourth-order Runge-Kutta approach. At the time level n + 1, the solution q n+1 
can be obtained as 


q <o > 

, (i> 

, (2 > 


n 

q 


, (0) 

q <°> 


- 4t 
3J 


At 

6J 


(Fq 


( 1 ) 




) 

) 
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( 14 ) 


( 4 ) 

q 

n+l 

q 


< (0) - £ <F,«> - „<»>) 
-(4) 


where At is the time step. 


2. 1 • 5 Boundary Conditions 

The boundary conditions on the solid boundaries are 
v • n = 0 

In other words, there is no convected flow across the solid surfaces. The 
far-field boundaries are posed as free-stream conditions with one-dimensional 
gas equation correction. 


2.1.6 Solution Iterations 

The flow properties are initialized as 


/ \ 

/ \ 

p 


o 

Q. 

pu 


p u cosa 

r O a> 

pv 

* 

P 0 u » sina 

pw 


0 

oe 


p e 


k ° ° 


where a is the angle of attack. For time-dependent calculation, the 
properties can be updated after each finite volume and Runge-Kutta 
iteration. For steady-state calculation, no time differencing Is needed. 
Therefore, only the finite-volume Iteration is used. 

2.1.7 Thrust Vectoring Application 

For a thin jet stream, the relation between pressure difference and the 
jet-path curvature can be expressed as 
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( 16 ) 



r 3r 


This is obtained by balancing the centrifugal force with the pressure force. 
Solving for the radius of curvature, it can be obtained as 


Ap 


(17) 


whe 


as 


re tj is the jet stream thickness. Also, the curvature 1/r can be written 


r 


li 


[l + ( Z :) 2 1 3/2 


2 


p jYj 


(18) 


The thrust vectoring application with FLO-59 is based on Equation (18). 
Initially, the calculation starts with an assumed jet shape. The grid has to 
be constructed following this initially guessed path. Then from the flow 
solver, pressure difference across the jet can be obtained. A new jet path 
can be updated with the relation described in Equation (18). Another set of 
grids should be built based on the new jet path. The iteration goes on until 
the jet path does not change. 

This approach is suitable for thin-jet calculation. For a thick jet, the 
formulation of Equation (17) is not valid. Another drawback of this approach 
is that the nozzle exit geometry cannot be accurately modeled. The main 
reason for this drawback is the limitation of single-block grids. The geo- 
metric unwrapping operation in FLO-59 is along the streamwise wing surface. 
Following this orientation, there would be a geometric discontinuity at the 
nozzle exit. This problem can be solved with the multiblock operation. 
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2.2 CFL3D (Thin-Layer Navier-Stokes Method) 


2. 2. 1 Governing Equations 

From conservation of mass, momentum, and energy, the Navier-Stokes 
equations can be written as 

I? + fc (F ' F - ) + (G - G J + fer (H - H_) - 0 


where 

/ P ' 
pu 
pv 
pw 

^ / 

Flux in x-direction: 


Q = 


F = 


by 


b z 
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N, 

pu 
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pu 

+ p 


T xx 

puv 


F v - 

T xy 

puw 



T xz 

(e + 

p)u 


UT + VT + WT - q 

X 

J 


xx xy xz 


( 19 ) 


( 20 ) 


Flux in g-direction: 

\ 
pu 

puv 
2 . 

pu + p 
puw 

(e + p)v 






G v = 
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T yx 

T yy 
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UT 
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+ VT + wr - q 
yx yy yz 4 y 




(21) 
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Flux in z-direction: 


/ 

\ 


i' n 

pw 



0 

pwu 



X 




zx 

pwv 


H - 

X 



V 

zy 

2 




pw 

+ p 


X 




zz 

(e + 

p)w 


UT + VT + WT - q 

\ 

/ 


zx zy zz 


( 22 ) 


2.2.2 Coordinate Transformation 


Because of the irregular grid spacing in the physical domain (x-, y-, z- 
coordinates), the order of computational accuracy is difficult to predict. 
From a scheme-stability point of view, the computation will be better per- 
formed in a generalized coordinate domain (so-called computational domain). 
The generalized coordinate transformation can be described as 


£ * S(x, y, z, t) 


T) » ri(x, y, z, t) 


C = C(x, y, z, t) 

(23) 


The Jacobian of this generalized transformation can be obtained as 

j = 3(S> T1, Q 

d(x, y, z) 

' We ■ y c z n > " VVt ' x eV * Wc ' Vc >r ‘ 

(24) 

The governing equations become 


a 

22 

dt 


> A ^ A A AAA. 

+ f 7 (F-F)+f-(G-G)+f— (H-H) = 0 

0£ OT) V 5t V 


(25) 
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where 


~ * I yr I A a a 

F - F v ■ J R* (F - F v )+ V® ' V + V” • V + «t Ql 

G - G „ - h x < F - f v ) + i y w - <=„> * i z CH - C v ) + r) t Q] 

* A [yr I A * a 

H - H v ■ J K, (F - V + Cy « - v + C z « - V + c t Q] 

<V K y , K 2 , K t ) - (K x , K y , K z K t )/|VK| 

|7K| - (K 2 + K 2 + K 2 ) 1/2 

1 x y C 

K 3 (5» T) » C), respectively, (26) 


2.2.3 Finite-Volume Scheme 


The governing equations, while written in generalized coordinates, are 
used in the finite volume formulation. Integrating equations over control 
volume bounded by lines of constant in Tp, C"directions yields 


A A 


A A 


+ (F ‘ Vi+l/2,j,k (F F v ) i-l/2,j,k 


A A 


+ (G * G vh , j+l/2,k - < G - ®v>l,]-l/2,k 


A A 


+ «- Vi.j.k+1/1 - (H - Vu.k-uz -° 


( 27 ) 


2.2.4 Thin-Layer Navier-Stokes Approximation 

It is known that the velocity and pressure vary rapidly in the direction 
normal to the surface. On the other hand, flow characteristics are generally 
smooth along the surface, unless singularities exist. Retaining only the 
viscous terms corresponding to the direction normal to the surface can save 
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computing resources, yet maintaining enough computational accuracy. In this 
manner, the thin-layer Navier-Stokes equations can be written as 

if + fe <?) + h (q) + tc ( " ' V - 0 < 28 > 

Certainly, viscous terms can also be retained in l , q directions because, in 
CFL3D, solid surfaces can be placed on any block surfaces. The laminar vis- 
cous terms can be included in all 5, q, C directions simultaneously. However, 
turbulence viscous terms can only be added at most in two directions simul- 
taneously. 


2.2.5 Approximate Factorization 


The thin-layer Navier-Stokes equations can be solved by approximate 
factorizations. CFL3D has the options of two— factor and three— factor approx- 
imate factorizations. Taking the three— factor approximate factorization, for 
example, at time step n + 1, Equation (28) can be written as 


r 1 + s &£ . «• <12. , r /d H _ 9Hv M n An D n 
JAt 3Q 6 n 3Q + VdQ SQ )] AQ = R 

where R n is the residual. Equation (29) can be factorized into three 
equations along £ , q , C sweep directions. 


(29) 


[ 

[• 

[■ 


I 


JAt 

I 

JAt 

I 

JAt 


. 3F,n * n 
+ gp] AQ - -R 


, . aG,n ** , I * 

+ \ 8Q 1 4Q ’ 


3H 

+ 5 c ( ^-^ )]nAQ = ( k )AQ 


(30) 

(31) 

(32) 
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The solution can be updated as 

Q n+1 * Q n + AQ (33) 

2.2.6 Upwind Differencing Scheae 

Two distinct methods are provided in this program, i.e. flux difference 
splitting and flux vector splitting. 


2.2.6. 1 Flux Difference Splitting: The generalized fluxes F, G, and H are 

split into forward and backward contributions according to the signs of eigen- 
values of the Jacobian matrices. For example, considering the flux in the £ 
direction, 


A - = TAT' 1 - T(A + + A )T 1 

o Q 


(34) 


where A 


± _ \ ± lxl 


, A are diagonal matrices formed from the eigenvalues of A; 


/v 

i.e., A = diag( U ^ 


u|Vg 

J 


JV /N» 

u|ASl (u + a)|vgl (u - a)!vgK 
J J ’ J ' 


and the symbol indicates a Roe-averaged variable. The left-hand side of 
Equation (30) can be expressed as 


'3ST + 6 ; H'-m* ■ + V + + « 5 ) 

In the same manner, the flux in the tj- and C-directions can be split 
accordingly. 
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2*2* 6* 2 Flux Vector Splitting: The generalized fluxes F, G, and H are split 

into forward and backward contributions according to the flow speed* 

For example, considering the flux in the £ -direction in Equation (28): 


A | i A __ 

6^F ■ 5^F + 6^F 


(36) 


where 6^ and 6^ denote backward and forward difference. 
Define : 


■ u/a and u * u/|V£| 


(37) 


for supersonic flow (|M^| > 1) 


F = F , F =0, M^>1 


A A 


F * 0, F - F, M < - 1 


(38) 


for subsonic flow (|M^| < 1) 


LZSj 


mass 


f [k (-u ± 2a)/y + u] 
mass x T 


f [k (-u ± 2a)/y + v] 
mass y 1 


f [k (-u ± 2a)/v + w] 
mass z 


energy 


mass 


±pa(M^ 


± 1)^/4 


<M ? ' a> 


f energy " f Lss ( [ " <Y ' 1);2 1 2<Y ' ^ + 2 ^ ] ° + ^ + ^ + ” 2)/2 > 


, k y - y|V5| , k, - « t /|Tt| 


(39) 


The flux in the ri - and £ “directions can be obtained in the same manner. 
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2.2.7 Multiblock Operation 


For complicated configurations, the traditional single-block flow solver 
cannot predict the flow properties without sacrificing accuracy at some por- 
tion of the flow field. For example, the unwrapping type grid system can only 
describe a simple configuration, such as wing alone or body alone cases. 
Applying the same grid system to a wing-body configuration, grid smoothing is 
needed. Also, this grid system usually produces coarse grids on the fuselage 
while using the same number of grid points on the wing in the unwrapping 
direction. For a configuration with integrated nozzles, the single-block 
calculation (grid generation and flow solver) cannot produce accurate re- 
sults. Therefore, raultiblock calculation becomes inevitable in this research. 

In the multiblock calculation, the computational domain is divided into 
several hexahedral blocks. Basically, the program solves the thin-layer 
Navier-Stokes equations in each individual block. Boundary conditions on each 
boundary surface of each block has to be posted. If a boundary surface Is 
designated as a patch surface, flow properties will be provided from neighbor- 
ing blocks connected to this patch surface. With this useful technique, 
complicated geometry can be solved without smoothing out the geometry or 
losing accuracy. 

2.2.8 Jet Application 

In the program CFL3D, eight standard boundary conditions are available. 

A list of the standard boundary conditions can be found in the Appendix. 
Special purpose boundary conditions can be added following the multiblock 
coordinate orientation. In the current work, a jet block has to be formed to 
simulate jet exit conditions. The jet exit plane is treated as one boundary 
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surface of the jet block* A special routine has been developed to specify the 
jet exit conditions and include the jet in the computation. On the other 
boundary surfaces, boundary conditions can be specified as patching or in/out 
flow surfaces. Examples of designating boundary conditions on block surfaces 
will be included in this report. 

For thrust vectoring applications, to separate the jet momentum from the 
free stream, jet blocks should be constructed following the jet trajectory. 

The jet path depends on the strength of the jet as well as the free-stream 
velocity. Investigation of the jet deformation shape and effects of jet 
conforming grids will be continued In the next period of this study. 

2. 3 Ffrmerical Grid Generation 

Generally, flow equations describe the flow physical properties; and the 
grids represent the discretized flow field that the flow equations influence. 
Transformation type grid generation can be used to describe simple geometries. 
However, for complicated configurations, numerical grid generation becomes an 
indispensable tool. Currently, the program EAGLE (Eglin Arbitrary ^Geometry 
Implicit Euler) has been used for this thrust vectoring research. Features of 
this program Include two-dimensional surface generation and three-dimensional 
grid generation. In CFD applications, an Individual computational domain can 
be seen as a hexahedral block. The surface generation produces coordinates on 
each boundary surface and/or any surface in the block. Consequently, the 
three-dimensional grid generation collects information from two-dimensional 
surface generation and produces the entire grids. Two grid generation mecha- 
nisms, algebraic and elliptic generation systems, are available in this program. 
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2*3*1 Surface Generation 


If the boundaries of a surface are specified, the coordinates on the 
surface can be obtained by various interpolation methods. Popular interpo- 
lation methods include Lagrange interpolation, Hermite interpolation, spline 
Interpolation, etc. In the EAGLE program, the Lagrange interpolation is 
recommended, although other interpolation functions are available. 

The general form of one-dimensional (£) Lagrange interpolation is 


i,- ! 5 


r<5) - i ♦„<!*„ • ♦„<?> ■ , n , 

n a l 1=1 ^n 


(40) 


where N is the number of point in the £ interpolation domain and 
£ e [0, I] . For example, for a two-point Lagrange interpolation, 


Y(0 - (1 - \)y + | y 2 


(41) 


The general form of two-dimensional (£, r\) Lagrange interpolation is 


Y(S, tj) = I «D n (f)Y(C n , ti) + l ti m ) 

n=l m=l 

N M 

- I I » T) ) 

, , nimJ n m 

n=l m=l 

where c|> (*^r) = II — 

" J l-l "m ' "t 


(42) 


M is the number of points in the -p interpolation domain, and 
T) e [0, J]. 

From Equation (3), the coordinates of the grid points can be obtained. 
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2*3*2 Three-Dimensional Grid Generation 


2* 3* 2.1 Algebraic Grid Generation: In a three-dimensional domain (£ , t ) , C)> 

the Interpolation becomes 


Y(5, t), C) 


N 

I 

n=l 




M 


"• C) + Ij I.- 5) 


m ; 




The grids can be produced by using Equation (43). These are the so- 
called algebraic grids. Usually this type of grid generation system produces 
smooth grid coordinates. The main advantages of this algebraic grid gener- 
ation are simplicity and saving of computing time. From Equation (43) the 
grids are produced based on the boundary spacing. Therefore, boundary point 
spacing is critical in this system. For example, in flow calculation the 
spacing near abruptly changing contour should be small. Therefore, the 
spacing should be small while specifying the boundary points for surface 
generation in EAGLE. In other words, users need good judgement while using 
the algebraic grid generation system. One way to avoid the disadvantages 
inherent in using the algebraic grid generation system is to use the elliptic 
grid generation system. 


2. 3.2.2 Elliptic Grid Generation: The purposes of elliptic grid generation 
system are controlling coordinate line distribution, and orthogonality in the 
field. Elliptic grids can be generated by solving the Poisson equations: 
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vV - P 1 (i - 1, 2, 3) (44) 

in which the control function is to control the spacing and orthogo- 
nality. Detailed formulation of the control functions can be found in Ref- 
erence 7. While solving the Poisson equations, the algebraic grids obtained 
from transfinite interpolation are used for initial guesses. Then an iter- 
ative procedure is used to solve Equation (44) until satisfactory results are 
obtained. 
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3* RESULTS 


Five configurations have been tested in this study. The first three test 
cases are configured to test the FLO-59 program. They are 

(1) A rectangular wing integrated with a nozzle 

(2) A 30-degree swept wing integrated with a nozzle 

(3) A 58-degree delta wing. 

To verify the new routines added in CFL3D, two rectangular nozzle configu- 
rations are tested. 

3» 1 Rectangular Wing Integrated with a Nozzle 

As shown in Figure 1, this rectangular ONERA wing (Ref. 8) is modified 
with the jet nozzle distributed between 40% and 70% of the semispan. The 

thickness of the nozzle is 4% of the root chord. The jet is deflected at an 

\ 

angle of 30 degrees relative to the chord. The thrust coefficient is 0.2. 
Pressure distribution at 25%, 45%, and 75% of the semispan is shown in Figure 2. 

3.2 30-Degree Swept Wing Integrated with a Nozzle 

This configuration also has the ONERA wing section distribution (Fig. 

3). The jet nozzle is distributed between 40% and 70% of the semispan. The 
thickness of the nozzle is 4% of the root chord. The jet is deflected at an 
angle of 30 degrees relative to the chord. The thrust coefficient is 0.2. 

Two different grid systems are tested with the same jet exit condition. The 
first grid system is constructed with the cells following the jet path; i.e., 
the cut line of the mesh being identical to the jet path centerline. In this 
manner, the calculation direction follows the physical flow. Hence, the Kutta 
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condition can be satisfied. For the purpose of comparison, another grid 
system with a straight cut line representing a straight wake has also been 
tested. The pressure coefficient distributions are shown in Figure 4. It is 
seen from the calculated results that the jet effect is diminished in the 
straight-wake case. At present, no experimental data for this test case are 
available. However, the results from the deformed wake formulation represents 
a more realistic solution. Therefore, it is concluded that the grids should 
follow the jet path to correctly model the jet effect. 

3.3 A 58-Degree Delta Wing 

One of the objectives in this study is modeling a wing-body configuration 
with a vectored jet, as shown in Figure 5. This aerospace plane has a highly 
swept lifting surface for high-speed flight conditions. Highly swept wings 
have the tendency to produce vortices. As known, the vortex is mainly gener- 
ated by pressure difference. In the situation of thrust vectoring, the pres- 
sure difference across the afterbody is largely increased. Therefore, even at 
low angle of attacks, vortex flow is likely to occur. 

A 58-degree swept delta wing is modeled to investigate the leading-edge 
flow (Fig. 6). Pressure distribution at 25%, 45%, and 75% are shown in Figure 
7. The leading-edge vortex is sensitive to the grid spacing near the leading 
edge (Ref. 9). In the calculation, the grid spacing is very dense near the 
leading edge. Reference 9 also demonstrates the success of capturing leading- 
edge vortex with the program CFL3D. 
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3. A A Rectangular Nozzle without Turning Vanes 


A new routine has been added in the program CFL3D to simulate the jet 
thrust. The jet exit plane is modeled as a block boundary surface. The jet 
velocity, temperature, and pressure can all be specified on this plane, A 
rectangular nozzle without turning vanes is modeled to test this new routine. 
This nozzle is featured as a 2D nozzle with an aspect ratio of 4,6, As 
shown in Figure 8, the jet exit Mach number is 0,9 with the free-stream Mach 
number equal to 0.1. No similar experimental data are available for compar- 
ison. However, from the theory of jet expansion, the calculated results match 
the entrainment and expansion concept. In this case, only the jet velocities 
are specified at the nozzle exit. For comparison with a real jet, the jet 
pressure as well as energy should also be specified. 

3*5 A Rectangular Nozzle with Thrust Vectorin g Effect 

\ 

This configuration is designed to test the thrust vectoring concept. As 
shown in Figure 9, the nozzle exit is deflected at 30 degrees. Figure 9 shows 
the block arrangement and specified boundary conditions on the block surfaces. 
A jet block is formed following the jet trajectory, which contains the jet 
momentum. From multiblock formulation point of view, the boundary surface 
positions of this jet block are not critical because the interpolations among 
blocks can capture the momentum and energy transition. However, the accuracy 
of the interpolations critically depends on the grid orthogonality and spacing 
on both sides of the interface. To eliminate the possibility of interpolation 
error, the jet block is constructed containing most of the jet momentum in the 
region near the exit. 
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The jet exit Mach number is 0.8, and the free-stream Mach number equals 
0.25. The angle of attack in this situation is zero. Again, no similar 
experimental data are available for this configuration. The symmetric plane 
u-w velocity profile shown in Figure 10 does demonstrate all features of a 
deflected jet. Also, pressure distribution in the longitudinal direction is 
shown in Figure 11. The sectional lift coefficient due to the deflected 
thrust is approximately 1.0. Under the same situation, the lift coefficient 
should be nearly zero without thrust vectoring effect. The calculation demon- 
strates the concept of induced lift force due to the deflected thrust. Owing 
to the pressure difference across the jet stream, vortices are expected to 
occur. Figure 12 shows the v-w velocity profile in a cross flow plane. The 
cross-flow plane is located at ten nozzle heights from the exit. As expected, 
vortices do occur in the jet stream. 
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4. CONCLUDING REMARKS 


The programs FLO— 59 and CFL3D were used to predict the thrust vectoring 
effects. Computational grids are generated by the program EAGLE, Complicated 
geometries, such as nozzle and conforming jet boundaries, can be modeled 
accurately by this program. 

The Euler flow solver FLO-59 produces reasonable results for wings inte- 
grated with jet-flap type exit. The pressure difference due to thrust vector- 
ing near the trailing edge can be reasonably calculated by this program. For 
afterbody nozzle configurations, the raultiblock flow solver CFL3D must be 
used. The raultiblock function enables users to model the nozzle in detail. 
Computational results demonstrate the flow field asociated with a deflected 
jet. Also, the induced lift can be predicted by CFL3D. 

The computational results presented in this report demonstrate the capa- 
bility of the flow solvers and the grid generation program. These tools will 
be applied to the study of wing-body configurations with thrust vectoring in 
the remaining period of this study. 
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Figure 1: 


Grids on the Rectangular Wing and Jet-Path Surface. 
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Pressure Distribution 



(o) 2y/b = 0.25 


Figure 2: Rectangular Wing Chordwise Pressure Distribution at M 
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(c) 2y/b = 0.75 


Figure 4. Continued 
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Pressure Distribution 



(d) 2y/b = 0.25 
Figure 4. Continued 
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Figure 4. Concluded 
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Pressure Distribution 


Figure 7 
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(b) 2y/b = 0.45 
Figure 7: Continued 
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Pressure Distribution 



(c) 2y/b = 0.75 
Figure 7: Concluded. 
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Figure 8: 


Section Velocity Profile of a Straight Jet. 

M - 0.1; M - 0.9; 2y/b - 0. 

» jet 
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Figure 10: Sectional u-w Velocity Profile of a Jet Deflected at 30 Degrees. 

2y/b - 0; - 0.25; M Jet - 0.8. 
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Pressure Coefficient 
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Figure 11: Pressure Distribution along a Thrust Vectoring Nozzle. 

2y/b - 0; - 0.25; M jet * 0.8. 
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Figure 12: Cross Plane v-w Velocity Profile at x /^ nozz i e “ 10 from the Exit. 


47 




APPENDIX 


Standard Boundary Conditions: 


1000 

free stream 

1001 

reflection x-z plane 

1002 

extrapolation 

1003 

inflow/outflow 

1004 

viscous surface 

1005 

inviscid surface 

1006 

x-y plane 

1007 

axi symmetric 


Special Purpose Conditions: 
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jet exit plane 


